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We report on an all-sky search with the LIGO detectors for periodic gravitational waves in the 
frequency range 50-1100 Hz and with the frequency's time derivative in the range —5 x 10~ 9 - 
Hz s _1 . Data from the first eight months of the fifth LIGO science run (S5) have been used 
in this search, which is based on a semi-coherent method (PowerFlux) of summing strain power. 
Observing no evidence of periodic gravitational radiation, we report 95% confidence-level upper 
limits on radiation emitted by any unknown isolated rotating neutron stars within the search range. 
Strain limits below 10 -24 are obtained over a 200-Hz band, and the sensitivity improvement over 
previous searches increases the spatial volume sampled by an average factor of about 100 over the 
entire search band. For a neutron star with nominal equatorial ellipticity of 10 -6 , the search is 
sensitive to distances as great as 500 pc — a range that could encompass many undiscovered neutron 
stars, albeit only a tiny fraction of which would likely be rotating fast enough to be accessible to 
LIGO. This ellipticity is at the upper range thought to be sustainable by conventional neutron stars 
and well below the maximum sustainable by a strange quark star. 

PACS numbers: 04.80.Nn, 95.55.Ym, 97.60.Gb, 07.05.Kf 



I. INTRODUCTION 

We have carried out an all-sky search with the LIGO 
(Laser Interferometer Gravitational-wave Observatory) 
detectors jU |3] for periodic gravitational waves, using 
data from the first eight months of LIGO's fifth science 
run (S5). We have searched over the frequency range 
50-1100 Hz, allowing for a frequency time derivative in 
the range —5 x 10~ 9 -0 Hz s _1 . Rotating neutron stars 
in our galaxy are the prime target. At signal frequencies 
near 100 Hz we obtain strain sensitivities below 10 -24 , 
a strain at which one might optimistically expect to see 
the strongest signal from a previously unknown neutron 
star according to a generic argument originally made by 
Blandford (unpublished), and extended in our previous 
search for such objects in S2 data [3]. A recent refine- 
ment of the argument jl] gives less optimistic estimates, 
but these too are surpassed by the experimental results 
presented here. 

Using data from earlier science runs, the LIGO Scien- 
tific Collaboration (LSC) has previously reported on all- 
sky searches for unknown rotating neutron stars (hence- 
forth designated as "pulsars" here). These searches have 
been performed using a short-period coherent search in 
the 160.0 - 728.8 Hz frequency range [3] , and using a long- 
period semi-coherent search in the 200 - 400 Hz frequency 
range in the S2 data [5| and the 50- 1000 Hz range in the 
S4 data [6]. Einstein@Home, a distributed home comput- 
ing effort [7], has also been running searches using a co- 
herent first stage, followed by a simple coincidence stage, 
for which S3 and S4 results have been released [3 13] ■ 

The data collected in the S5 data run were more sensi- 
tive than in previous data runs, and the amount of data 
used here is an increase by a factor of eight over that 
reported from the S4 data run [6], resulting in upper 
limits on periodic gravitational waves about a factor of 
3-6 lower than those from the S4 data, depending on 
source frequency. This improvement gives an increase 
in sampled galactic volume by about a factor of 100, de- 
pending on the assumed source frequency and spin-down. 



At a signal frequency of 1100 Hz we achieve sensitivity 
to neutron stars of equatorial ellipticity e ~ 10 -6 at dis- 
tances up to 500 pc (see [6 j for relations). This ellipticity 
is at the upper range thought to be sustainable by con- 
ventional neutron stars [10 and well below the maximum 
sustainable (10 -4 ) by a strange quark star [11] . The num- 
ber of undiscovered, electromagnetically quiet neutron 
stars within 500 pc can be estimated to be O(10 4 — 10 5 ) 
from the neutron star birth rate [IS], although it is likely 
that only a tiny fraction would both be rotating fast 
enough to be accessible to LIGO [T3] and remain in the 
local volume over the age of the galaxy [2] . Only ~ 25 
radio or x-ray pulsars have been discovered so far within 
that volume [i~5] . 



II. THE LIGO DETECTORS AND THE S5 
SCIENCE RUN 

The LIGO detector network consists of a 4-km interfer- 
ometer in Livingston Louisiana, (LI), and two interfer- 
ometers in Hanford Washington, one 4-km and the other 
2-km (HI and H2). 

The data analyzed in this paper were produced in 
the first eight months of LIGO's fifth science run (S5). 
This run started at 16:00 UTC on November 4, 2005 
at the LIGO Hanford Observatory and at 16:00 UTC 
on November 14, 2005 at the LIGO Livingston Observa- 
tory; the run ended at 00:00 UTC on October 1, 2007. 
During this run, all three LIGO detectors had displace- 
ment spectral amplitudes very near their design goals of 
1.1 x 10 _19 mHz -1 ^ 2 [16 in their most sensitive frequency 
band near 150 Hz. (In terms of gravitational- wave strain, 
the H2 interferometer was roughly a factor of two less 
sensitive than the other two; its data were not used in 
this search.) 

The data were acquired and digitized at a rate of 
16384 Hz. Data acquisition was periodically interrupted 
by disturbances such as seismic transients (natural or 
anthropogenic) , reducing the net running time of the in- 
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terferometers. In addition, there were 1-2 week commis- 
sioning breaks to repair equipment and address newly 
identified noise sources. The resulting duty factors for 
the interferometers were approximately 69% for HI, 77% 
for H2, and 57% for LI during the first eight months. A 
nearby construction project degraded the LI duty factor 
significantly during this early period of the S5 run. By 
the end of the S5 run, the cumulative duty factors had 
improved to 78% for HI, 79% for H2, and 66% for LI. 
For this search, approximately 4077 hours of HI data and 
3070 hours of LI data were used, where each data seg- 
ment used was required to contain at least 30 minutes of 
continuous interferometer operation. 



III. SIGNAL WAVEFORMS 

The general form of a gravitational-wave signal is de- 
scribed in terms of two orthogonal transverse polar- 
izations defined to be "+" with waveform h+(t) and 
"x" with waveform h x (t), for which separate and time- 
dependent antenna pattern factors F + and F x apply, 
which depend on a polarization angle ip [17] . For peri- 
odic gravitational waves, which in general are elliptically 
polarized, the individual components h +tX have the form 
h + (t) = A + cos$(t) and h x (t) = A x sin$(t), where A + 
and A x are the amplitudes of the two polarizations, and 
$(t) is the phase of the signal at the detector. For the 
semi-coherent method used in this search, only the in- 
stantaneous signal frequency in the detector reference 
frame, 2wf(t) = d$(t)/dt, needs to be calculated. For an 
isolated, precession- free, rigidly rotating neutron star the 
quadrupolar amplitudes A + and A x are related to wave 

amplitude, h , by A + = h +c ° s - and A x = ho cos i, 
where i describes the inclination angle of the star's spin 
axis with respect to the line of sight. For such a star, the 
signal wave frequency, /, is twice the rotation frequency, 

fr- 

The detector reference frame frequency f(t) can, to 
a very good approximation, be related to the frequency 
f(t) in the Solar System Barycenter (SSB) frame by [5] 
f(t) - f(t) ~ f{t)^^, where v(i) is the detector's ve- 
locity with respect to the SSB frame, and n is the unit- 
vector pointing from the detector toward the sky location 
of the source [5] . 



IV. ANALYSIS METHOD 

The PowcrFlux method used in this analysis is de- 
scribed in detail elsewhere [6] and is a variation upon 
the StackSlide method [18]. Here we summarize briefly 
its main features. 

A strain power estimator is derived from summing 
measures of strain power from many short, 50%-overlap, 
Hann-windowed Fourier transforms (SFTs) that have 
been created from 30-minute intervals of calibrated strain 



data. In searching a narrow frequency range (0.5 mHz 
spacing) for an assumed source sky location, explicit cor- 
rections are made for Doppler modulations of the appar- 
ent source frequency. These modulations are due to the 
Earth's rotation and its orbital motion around the SSB, 
and the frequency's time derivative, /, intrinsic to the 
source. Corrections are also applied for antenna pattern 
modulation, assuming five different polarizations: four 
linear polarizations separated by 7r/8 in polarization an- 
gle, and circular polarization. When summing, the vari- 
ability of the noise is taken into account with an SFT- 
dependent weight proportional to the expected inverse 
variance of the background noise power (see [5J [H] for 
detailed formulae). 

The search range for initial frequency /o values is 50 - 
1100 Hz with a uniform grid spacing equal to the size 
of an SFT frequency bin [1/(30 min)]. The range of / 
values searched is —5 x 10 _9 ~0 Hz s _1 with a spacing of 
5 x 10~ 10 Hz s _1 , since isolated rotating neutron stars 
are generally expected to spin down with time. As dis- 
cussed in our previous reports [3j [6] , the number of sky 
points that must be searched grows quadratically with 
the frequency / , ranging here from about five thousand 
at 50 Hz to about 2.4 million at 1100 Hz. The sky grid 
used here is isotropic and covers the entire sky. 

Upper limits calculated in this method are strict fre- 
qucntist limits on linear and circular polarization in small 
patches on the sky, with the limits quoted here being the 
highest limits in each 0.25-Hz band over broad regions 
of the sky. These are interpreted as limits on worst-case 
(linear polarization) and best-case (circular polarization) 
orientations of rotating neutron stars. Since the eight 
months of data analyzed here cover a large span of the 
Earth's orbit, providing substantial Doppler modulation 
of source frequency, contamination from stationary in- 
strumental lines is much reduced from earlier and shorter 
data runs. A total of only 0.6% of the search volume in 
sky location and spindown had to be excluded from the 
upper-limit analysis because of Doppler stationarity. 

The primary changes in the PowerFlux algorithm used 
in this search concern followup of outlier candidates. 
(The general method for setting upper limits is identical 
to that used in the S4 search (BJ.) Here we summarize 
the followup method used. Single-interferometer searches 
are carried out separately for the HI and LI interferome- 
ters, leading to the upper limits on strain shown in Fig. [I] 
and discussed below. During determination of the max- 
imum upper limit per sky region, per frequency band 
and per spin-down step, a "domain map" is constructed 
of local signal-to-noise ratio (SNR) maxima, with the do- 
mains ordered by maximum gridpoint SNR and clustered 
if close in direction and frequency. The 1000 domains 
with the highest maximum SNR are then re-analyzed 
to obtain improved estimates of the associated candi- 
date parameters, using a modified gradient search with a 
matched filter to maximize SNR with respect to source 
frequency, spin-down, sky location, polarization angle ip, 
and inclination angle t [19] . This maximization step sam- 
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pies frequency and spin-down much more finely than in 
the initial search. 

When all sky regions and all spin-downs have been 
searched for a given 0.25 Hz band for both HI and 
LI, the search pipeline outputs are compared, and the 
following criteria are used to define candidates for fol- 
lowup analysis. The HI and LI candidates must each 
have an SNR value greater than 6.25, and they must 
agree in frequency to within 1/180 Hz = 5.56 mHz, in 
spin-down to within 4xl0~ 10 Hz s _1 , and in sky loca- 
tion to within 0.14 radians. These conservative choices 
have been guided by simulated single-interferometer pul- 
sar injections. Coincidence candidates within 0.1 Hz of 
one another are grouped together, since most candidates 
arise from detector spectral artifacts that become appar- 
ent upon manual investigation. 

Candidates passing these criteria are subjected to a 
computationally intensive followup analysis that repro- 
duces the all-sky PowerFlux search in a 0.25 Hz band 
around each candidate, this time using the (incoherently) 
combined strain powers from both interferometers. Sky 
maps of strain and SNR are created and examined man- 
ually for each individual interferometer and for the com- 
bined interferometers. Spectral estimates from noise de- 
composition are also examined to identify possible arti- 
facts leading to the coincident outliers. 



V. RESULTS 

Figure [l] shows the lower of the HI and LI 95% 
confidence-level upper limits on pulsar gravitational wave 
amplitude ho for worst-case and best-case pulsar orienta- 
tions for different declination bands (each with different 
run-averaged antenna pattern sensitivity). As in the S4 
analysis, narrow bands around 60- Hz power mains har- 
monics, along with bands characterized by non-Gaussian 
noise, have been excluded from the displayed limits. Nu- 
merical values for frequencies and limits displayed in 
these figures can be obtained separately [20] . Systematic 
uncertainties on these values are dominated by calibra- 
tion uncertainty at the ~10% level. 

All outliers were checked for coincidence between HI 
and LI, as described above. In most cases single- 
interferometer spectral artifacts were readily found upon 
initial inspection, most of which had known instrumental 
or environmental sources, such as mechanical resonances 
( "violin modes" ) of the wires supporting interferometer 
mirrors, and power mains harmonics of 60 Hz. Other out- 
liers were tracked down to previously unknown electro- 
magnetic disturbances. For six coincidence candidates, 
no instrumental spectral artifacts were apparent. Their 
0.1-Hz bands and favored spin-down values are listed in 
Table [T] along with the maximum SNR's observed in HI 
and LI data. 

None of these six candidates was confirmed, however, 
as a detection of a constant-amplitude, constant-spin- 
down periodic source of gravitational radiation. In each 



Frequency band Spin-down 



(Hz) 


(Hz s- 1 ) 




HI SNR LI SNR 


867.2 


-4.3 x 10" 


-9 


6.27 


6.30 


941.0 


-2.0 x 10 


-9 


6.50 


6.67 


967.8 


-1.5 x 10 


-9 


6.26 


6.33 


979.5 


-5.0 x 10" 


-9 


6.40 


6.29 


1058.6 


-5.0 x 10" 


■10 


6.83 


6.38 


1070.2 


-3.0 x 10" 


■10 


6.72 


6.99 



TABLE I: List of coincidence candidates for which no instru- 
mental spectral artifacts were observed. 



case, we found that the combined Hl-Ll SNR did not 
increase by more than 0.6 (0.4) units over the minimum 
(maximum) of the single-interferometer SNR's, with four 
candidates showing a decrease for combined SNR. To un- 
derstand the expectation for a true signal, we carried out 
a posteriori software signal injections, which indicated 
that combined SNR should typically show an increase 
over minimum SNR by more than 2.0 units for a single- 
interferometer SNR threshold of 6.25. Hence we conser- 
vatively veto all candidates with an SNR increase less 
than 1 unit. In addition, manual exploration of these 
candidates was carried out, using larger portions of the 
S5 run's data, to determine whether SNR increased with 
additional data, and with subsets of the the original 8- 
month data, to determine whether a transient astrophys- 
ical source could explain the candidate. None of these 
explorations proved fruitful. 

We also note that multi-interferometer injections indi- 
cate that for signal frequencies above 850 Hz, the coinci- 
dence requirements in frequency and sky location could 
be tightened by a factor of five to 1 mHz and by a factor 
of seven to 0.02 radians, respectively with only a slight 
reduction in efficiency for true signals. None of the can- 
didates in Table [T] satisfies these tighter criteria. 

In summary, we have set strict, all-sky frequcntist up- 
per limits on the strength of continuous-wave gravita- 
tional radiation of linear and circular polarization, cor- 
responding to least favorable and most favorable pulsar 
orientations, respectively. Followup analysis of coinci- 
dence candidates with SNR > 6.25 did not yield a de- 
tection. The limits on detected strain can be translated 
into limits on equatorial ellipticity as small as 10 -6 for 
unknown neutron stars as far away as 500 pc. This ellip- 
ticity is at the upper range thought to be sustainable by 
conventional neutron stars and well below the maximum 
sustainable (10~ 4 ) by a strange quark star. 
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FIG. 1: Minimum (HI or LI) upper limits (95% CL) on pul- 
sar gravitational wave amplitude ho for the equatorial (red), 
intermediate (green), and polar (blue) declination bands for 
best-case (lower curves) and worst-case (upper curves) pul- 
sar orientations. Shown are all the minimum limits for each 
of the 11 spin-down values from —5 x 10 _9 Hzs _1 to zero in 
steps of 5 x 10~ 10 Hzs _1 . 
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